Spin Order accompanying Loop-Current Order in Cuprates 
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The theory of the long range order of orbital current loops in the pseudogap phase is generalized 
to include the effects of spin-orbit scattering. It is shown by symmetry arguments as well as by 
microscopic calculation that a specific in-plane spin-order must necessarily accompany the loop- 
current order. The microscopic theory also allows an estimate of the magnitude of the ordered spin- 
moment. Exchange coupling between the generated spins further modifies the in-plane direction 
of the spin moments. The structure and form factor for the spin and orbital moments combined 
with the induced spin order is consistent with the direction of moments deduced from polarization 
analysis in the neutron scattering experiment. 
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INTRODUCTION 

In all underdoped cuprates thermodynamic, transport, 
and spectroscopic measurements reveal the formation of 
a new state of matter which is commonly referred to as 
the "pseudogap" state. Any credible theory of cuprates 
must specify the nature of this state. One of the the- 
ories proposed predicts that this phase breaks time re- 
versal through ordered current loops in o-cu-o plaquettes 
without breaking translational symmetry. [l], [3] . The 
fluctuations about this state are such that there is no 
specific-heat singularity at the transition Qj. This re- 
moves the major difficulty in regarding that the pseudo- 
gap state represents a broken symmetry. The unit cell for 
YBCO is shown in fig.fTJA.). The predicted loop-current 
order in the copper oxide planes has the pattern shown 
in fig.|T|3)for one of the four possible domains. Evidence 
for such a state was obtained from ARPES using circu- 
larly polarized light [H in BISCCO, and more directly by 
recent polarized neutron scattering diffraction in YBCO 

a. 

While the spatial symmetry of moments of fig.([Tj3) is 
borne out by the neutron experiments, the direction of 
the magnetic moments is not consistent with the pre- 
dictions. The orbital moments should be normal to the 
o-cu-o plaquettes. The plaquettes are not co-planar with 
the two-dimensional cu-planes due to the buckling of the 
planes (see figdJA)); for YBCO in which the neutron 
scattering experiments were done the nearest neighbor 
o-cu-o plaquettes make an angle of about 7° with respect 
to the Cu-planes. Therefore a tilting of only about 7° of 
the moments with respect to the normal to the Cu-planes 
is expected. However this angle has been deduced to be 
45 ±20° [|. 

The purpose of this paper is to resolve this matter. 
The basic physical point we draw on [3| is that spin orbit 
interaction can lead to spin ferromagnetism in states with 
orbital currents 0]. We present general symmetry argu- 
ments supporting this and calculate microscopically the 
nature of spin order in YBCO for states of the symmetry 
consistent with the observations Q. The magnitude of 
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FIG. 1: (A)Crystal structure of YBCO and (B) the Current 
pattern in the observed time reversal violating states 



the spin-moment is estimated to be only 10 — 20% of the 
orbital moment. But the neutron scattering intensity de- 
pends on the spatial distribution or the form-factor of the 
moments besides their magnitudes. For the [Oil] Bragg 
peak studied in experiments, we find that the spin form 
factor is significantly larger than the orbital current form 
factor because the latter are spread out more inside a 
unit-cell than the former. The existing experimental re- 
sults may thus reconciled with the theory but a definitive 
confirmation awaits the measurement of the form-factors 
in experiments. 

The direction of the ordered in plane spin-moments is 
affected also by the exchange interaction between the mo- 
ments. We present a rough order of magnitude for this 
effect which suggests that the spin order may be quite 
complicated. Reliable theoretical estimates on the actual 
spin-order are very difficult to make at this point. How- 
ever, there are some general features of the results which 
are expected to be robust. The details of the magnetic 
order suggested as possible here can only be resolved in 
experiments which have greater accuracy than the one 
performed to date. A companion to this paper contains 
the details of the experimental results which were pub- 
lished as a short report earlier as well as an analysis of 
the data applying the ideas in this paper |T(| • 
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SYMMETRY CONSIDERATIONS 

We can deduce through arguments using the symmetry 
of the crystal structure, and of the orbital order param- 
eter, Fig. [TJ that spin-orbit coupling must induce an in- 
plane order of the spins. First we present an argument by 
constructing a free-energy invariant which is a product of 
the spin-orbit tensor, the orbital order and the possible 
spin-order. We will show that the same conclusions can 
also be obtained from general arguments patterned after 
those due to Dzyaloshinskii 111 ] . These arguments allow 
one to deduce the symmetry of the spin order parameter 
but the physical basis of these general arguments requires 
a microscopic theory which also allows us to obtain the 
magnitude of the effective moments. 

The crystal structure of YBCO is shown in figQjA). 
Notice that the Oxygen atoms are displaced from the 
plane formed by the Copper atoms; the copper oxide 
plane is not flat but dimpled. Cu's do not sit at a point 
of Inversion. This leads to a linear coupling between the 
spin and orbital degrees of freedom. 

Consider first the symmetries of the crystal structure. 
Due to the buckling of the Oxygens, the crystal structure 
in figlHA) 

(i) breaks inversion symmetry through Copper and Oxy- 
gen 

(ii) but preserves reflection symmetry about x, y, x + y 
and x — y. 

The spin-orbit coupling tensor A, which couples linearly 
the spin and the momentum of the electrons, respects 
these symmetries. Let us now look at the symmetries 
of the current loop order. The order parameter Mo, 
corresponding to the domain in figQjB), 

(i) breaks inversion symmetry through Copper and Oxy- 
gen 

(ii) preserves the reflection symmetry about (x — y) 

(iii) breaks reflection symmetries about (x, y, x + y) 
Let Ms specify the distribution and direction of pos- 
sible spin-order which obeys the translational symmetry 
of the crystal and other symmetries so that it must ac- 
company the orbital order. An invariant term in the 
free-energy of the form 



f so = A^MSM§ + M 
X 



(1) 



must then exist. Here \ is the spin susceptibility for the 
order specified by M$. Given eqnQ] 



(2) 



will be realized. 

Consider the symmetries that need to be satisfied by 
Ms. The product of A and Mo preserves inversion 
through Copper and Oxygen, but breaks time-reversal 
and the one mentioned reflection. Hence Ms must be 



(i) odd under reflection about (x, y, x + y) 

(ii) even under reflection about (x — y) 

(iii) even under inversion about Copper and Oxygen 

A spin-order Mg consistent with these requirements is 
shown for one of the two cu-o bilayers in figlU Under in- 
version through Copper, the spins remain the same thus 
satisfying inversion (since spin is an axial vector). All re- 
flections other than x = y are broken. Since the oxygens 
layers are below the cu-layers in one of the bi-layers and 
above it in the other, the spin-orbit coupling has oppo- 
site sign in the two bilayers. It follows that the direction 
of the moments specified by Ms is opposite in the two 
bi-layers so that the net moment per unit-cell is zero. 

In the argument above (AMq) acts as a net mag- 
netic field on the spins and the order is stabilized by 
the quadratic term in the free energy. One can also give 
an argument, which we find a bit more abstract, follow- 
ing Dzialoshinskii-Moriya (DM) [HE! for M s . In DM, 
one asks whether an anti-symmetric interaction between 
magnetic moments Ma and Mb of the form 



D AB ■ {M A x M B ) 



(3) 



has the symmetries of the lattice. The direction of Dab 
is specified by the crystal symmetry in relation to the 
position A and B of and Mb- The gene ral condi- 
tions on Dab have been given by Moriya 12j. We take 
M^ and Mb to be the moments at the position of the 
centroid of the two triangles with the currents shown in 
Fig. QJ5. Due to the buckling of the planes, there is no 
center of inversion in the vector connecting these two mo- 
ments. Then Dab ^ 0. A mirror plane perpendicular to 
AB bisects AB. Then Dab must be parallel to the mirror 
plane. There exists also twofold rotation axis perpendic- 
ular to AB which passes through the mid-point of AB. 
Then Dab must be perpendicular to the two-fold axis. 
Thus Dab is along the c-axis and in the mirror plane 
specified by its normal x — y passing through cu. Given 
such a Dab, a tilt of and Mb so that they have a 
finite in-plane is mandated by the term ([3]). This is con- 
sistent with the direction of the spin-order deduced from 
the previous argument and shown in Fig. O 

An important point to emphasize here is that the sym- 
metry considerations do not specify the relative orienta- 
tion of the spin on the Copper and Oxygen atoms. The 
moment in the unit-cell has a component in the x + y 
direction but the relative direction of the spin- moment 
on Oxygen and on Copper spin are not specified. The 
microscopic calculation in the next section based purely 
on spin-orbit scattering provides the result that these 
moments are parallel. However, once a spin-moment is 
generated by the "effective field" provided by spin-orbit 
coupling and orbital order, one must consider also the 
exchange interaction between them. The actual spin ar- 
rangement depends on the relative magnitude of the ex- 
change interactions to the " effective fields" for spin order 
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where a and (3 are spin indices, the sum is over the po- 
sition of the Copper atoms, t p d and t pp are renormal- 
ized hopping which include the effects of the large on 
site repulsions, and V is the Coloumb repulsion between 
the charges on the Copper and Oxygen. The symmetry 
breaking is captured by considering the mean field de- 
coupling of the quartic term. An operator identity for 
the interaction terms is: 
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FIG. 2: Calculated spin order in YBCO for the domain shown 
in figflTB) in the absence of exchange interaction 
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FIG. 3: Schematic spin order in YBCO for the domain shown 
in figlUB) in the presence of exchange interaction. The spins 
need not be collinear and the angle between the net oxygen 
spins (summed over all oxygens in an unit cell) and the Cop- 
per atom is dictated by the relative strength of the exchange 
and ordering fields. 



due to orbital order and spin-orbit coupling, 
change interactions will be considered later. 



MICROSCOPIC MODEL 



The ex- 



ier d and the two 
31. The unit cell 



croscopic model on the basis of the Cop 
Oxygen p x , y orbitals in each unit cell 
is shown in figO The Hamiltonian for the system is 

H ^ ^ j^W^za (.Pi+xa Pi — xa "h Pi-\-ya Pi— yet) (4) 

ia 

+ i pp (p\ +xa — Pi-xSj (Pi+ya — Pi-ya) 



d L d *»pl +x pPi+xf3 



in. 



i-\-xa.(3 | 



+ d ia dia +Pi+ x pPi+xP) • (5) 



Here ji^+xa/3 is the current-tensor between the sites i and 
i + x. 

3i,i+x<xf} = « (4 a Pi+xl3 - Pl +X f3 d ia) (6) 

The mean field ansatz is 



3i,i± 



^5, 



a(3 



Rexp(i(p) (7) 



So an order parameter in terms of ordinary current (and 
not spin-current) is sought. Symmetry requires if = tt/2. 

Due to the buckling of the planes there is finite overlap 
between the Oxygen (p x ,Py) and the Copper (d xz ,d yz ) 
orbitals respectively. The spin orbit interaction on Cop- 
per couples the d x 2_ y 2 orbital with the d xz and d yz or- 
bitals but such matrix elements are off diagonal in spin. 
Thus processes are allowed wherein the electron can hop 
from the ground state orbital of the Copper to that of 
Oxygen and also flip its spin in the process. The effec- 
tive hamiltonian generated by such processes is of the 
form [H [3 : 



i,S 



X 



i,i-\-S 



i+S0 



+ C.C. 



(8) 



where S = (±x, ±y), A M±X = A y and A M±y = 
— Xqx. The coupling constant is given by Ao = 
|^d x 2_ J/ 2| L \d xz ) t xz /e a \, where L is the angular momen- 
tum operator, t xz is the hopping matrix element between 
the d xz and p x orbitals and eo is the energy difference be- 
tween d x 2_ y 2 and d xz orbitals. The total Hamiltonian for 



The theory of loop-current order is derived from a mi- the system is 



H = H, 



mf 



Hp 



(9) 



We have included the spin exchange term in the Hamil- 
tonian, H ex which we will discussed later. The Mean 
Field Hamiltonian Q is obtained from from eqn[4] with 
the mean field ansatz (eqnJSJ) made to decouple the quar- 
tic interaction. 
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We first determine the spin state (in the absence of the in the basis {d^ , dj^ , p x ^ , _P x ki > ^yki > Pykt } ^ = H m j 
exchange coupling discussed below) and see that it indeed +Hso is given by : 
reproduces the results from general symmetry grounds 
obtained above. Fourier transforming the Hamiltonian, 



H 



2i (tpdS x (k) + Rc x (k)) I + i\oc x <7y 2i (t pd s y (k) + Rc y (k)) I - i\ c x cr x 

2i (i p ds x (k) + Rc x (k)) I + i\ c x cjy 4f PJ) s x (k) s y (k) I 



2% (tpdSy (k) + Rc y (k)) I - i\ c x a x 



it pp s x (k) s y (k) I 
I 







i!Pi 



where s x>y (k) = sin (fc,;a/2, fc y a/2), c Xi3/ (k) = 
cos {k x a/2, kyCi/2), I is the identity matrix, and ex's are 
the pauli matrices. Consider first no spin-orbit coupling, 
i.e. Ao = 0. This mean field Hamiltonian lead to the 
Time-reversal breaking of the loop-current phase with 
order parameter R, but for A = it preserves spin ro- 
tational invariance. The minima of the band is shifted 
from the T point corresponding to the fact that the 
ground state breaks inversion and the reflection sym- 
metry —x + y. The particular direction of wave-vector 
picked out by the ground state depends on the choice of 
domain. 



SPIN ORDER 

We now estimate numerically the direction and mag- 
nitude of the spin-moment in the absence of exchange 
coupling. To do so we discretize the Brillouin zone and 
for each wave vector, k, we find the six eigenstates and 
corresponding eigenvalue of (fTU|) . Given the eigenstates 
we can compute the contribution to the spin moment at 
the Copper and Oxygen sites by taking the expectation 
value of their respective spin operators. We than sum 
the contributions from all states below the chemical po- 
tential. For t p d = 1, t pp = 0.4, R = 0.1, and A = 0.1, we 
find, for the domain shown in figfTj3, that the moment 
is distributed as shown in figO To understand the ori- 
gin of the spin moment we plot in fig[4]the energy of the 
two topmost bands which are near half filling. For the 
occupied states of the bands the corresponding spins on 
the Copper atoms are shown in figE] All other bands are 
fully filled and do not contribute to total spin. By sum- 
ming over occupied states we get the net spin and from 
fig [5] we see that there is a net spin along the — x + y 
direction. 

Notice, as required by our earlier symmetry argument, 
the spin order breaks the x = — y reflection symme- 
try but not the x = y reflection symmetry. The mag- 
netic moment due to the orbital current is estimated as 
~ 7(a 2 /8) where I is the current and a is in-plane lattice 
constant. The current density is related to the order pa- 



{£} 




kx 



FIG. 4: Energy for the two topmost bands, plotted in 
Grayscale, and their respective Fermi energy contour. Split- 
ting shown is due to the presence of spinorbit coupling. Notice 
that band (b) is shallower than band (a) and that the minima 
of the bands is shifted from the F point reflecting the broken 
inversion symmetry. 



rameter as (R/i p d)evFCiQ 3 where vf is the fermi velocity 
and ao is the typical size of the atomic orbital. Thus the 
net moment is ~ (R/ipd){TimvF / £fo>v>)^b, where m is 
the mass of the electron and €f is the Fermi energy. For 
the band structure of YBCO, {hmvF / 'cf^o) ~ ^(1)> wc 
estimate the orbital moment to be 0.1/us for the values of 
parameters chosen. The magnitude of the spin-moment 
for the same parameters is estimated from the calcula- 
tions represented in Fig. to be 0.01/is on the Oxygen 
and Q.02fj,B on the Copper. 
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{a} 
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{b} 



FIG. 5: Orientation of spins on the Copper atom for the two 
topmost bands for the occupied states. 



effective field discussed above competes with the the cou- 
pling J. In particular one would expect a long wavelength 
modulation of the net moment depending on the relative 
magnitude of the spin orbit interaction and the exchange 
coupling. The in-plane moment is found to be of order 
10~ 2 /ib from which we can estimate the net effective field 
(AM ) ~ (Ms) Ep = lCr 2 eV/^ B . The exchange cou- 
pling Jcu-o is of order leV, while J cu - cu is of order O.lev. 
The effective exchange field defined as J cu -o (Ms) is then 
also of order 10 _2 eV / \ib- In figf3]we schematically draw 
the effect of exchange interaction on the in plane spins. 
Symmetries dictate the existence of in plane moments 
while the orientation of their ordering depends on the 
relative strengths of the exchange and ordering fields. It 
should also be clear that the competition between these 
effects will in general change the translational symmetry 
of the spin-pattern. Depending on the relative exchange 
parameters and the spin-orbit coupling, the spin-pattern 
can be very complicated and in general incommensurate. 
However, since these small moments are daughters of the 
orbital order, the latter is expected to be modified only 
weakly. At this point, it is not worthwhile to speculate on 
the details of the spin-order of the small spin- moment of 
O(10~ 2 )/xs since the exchange energies can be estimated 
only very imprecisely in the metallic or pseudogap state. 
We must rely on the details of the magnetic structure to 
be obtained from the neutron experiments but for such 
small moments, this is no easy task. 



As discussed in the section on symmetry, including ex- 
change changes the relative orientation of the Copper and 
the net Oxygen spins and might also affect the magni- 
tude. 



Effect of Exchange Interaction on the Order of Spins 

The microscopic calculation above provides the result 
that the net spin on the two oxygens and on cu, generated 
due to spin-orbit coupling, are parallel. However, once 
a spin- moment is generated by the "effective field" pro- 
vided by spin-orbit coupling and orbital order, one must 
consider also the exchange interaction between them to 
determine the direction of the moments. Exchange intro- 
duces a coupling whose leading terms are of the form 



— Jcu-o Sj • (S° +x + S t 



i—yj 



J cu- 



fin 



<%,j> 



where J cu -o and J cu -cu are the exchange couplings be- 
tween nearest neighbor cu and oxygen and nearest neigh- 
bor cu's respectively. We expect the exchange interaction 
to be antiferromagnetic. Hence the induced order via the 



EXPERIMENTS 

From the symmetry and microscopic analysis above, 
we have shown that there are two sources of modulated 
magnetic fields within the sample: current loops and spin 
order. Using polarized neutron scattering, Fauque et 
al.[6|, performed a detailed study of five different sam- 
ples of YBCO (four underdoped and one overdoped) to 
look for magnetic ordering. The most detailed measure- 
ments are on an untwinned sample where the uncertainty 
is smaller. 

The principle conclusions reached were the following: 
(i) A new magnetic contribution to scattering intensity 
arose at the [Oil] Bragg peak in all underdoped samples 
below a temperature which increases as the sample is 
progressively underdoped. (ii) No signal was seen at the 
[002] Bragg peak. 

(iii) No new Bragg peaks appear ruling out breaking of 
translation invariance. 

(iv) The data could be fit to the current loop model pro- 
vided one assumed that the moments were located at the 
centroid of the triangles. 

(v) Assuming that all the signal was due to a single source 
of magnetic ordering, i.e. it has a unique form-factor 
and structure factors, the moment had to be tilted away 
from the c-axis. The angle was largest for the detwinned 
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sample with the moment being at ~ 45° ± 20°. Our 
finding that the in-plane moment is due to spins while the 
out of plane component is due orbital order necessitates 
a reevaluation of these numbers. 

To fit the neutron scattering data one has to assume a 
model for the magnetic moments in the system. From the 
observations above we conclude that within the experi- 
mental uncertainties, (~ 0.01^^), the moment is com- 
mensurate with the lattice and that the net moment in 
a copper oxide plane is zero. The latter follows from 
absence of observable magnetic signal at [002] . For loop- 
current order, the spin flip signal should appear at Bragg 
Peaks [0, K, L] , [H, 0, L], [H, H, L] and [H, -H,L]. The 
dimpling of the plane implies that these moments are at 
an angle of ~ 7° with the c-axis. To understand the 
origin of the larger deduced angle as stated in point (v) 
above we have to take into account the fact that the in- 
plane moments arise due to ordering of spins while the 
out of plane component is due to current loops. The cor- 
responding moments have very different form and struc- 
ture factors. For any given ordering of moments M(r) the 
spin flip scattering intensity at Bragg peak Q = [H, K, L] 
for polarization of the incident Neutron parallel to Q is 




/(q) ex |/ Q | 2 |S(Q)| 2 |A/J 2 (12) 
M ± = Q x (M x Q) /Q 2 

where /q is the form factor, given by the fourier trans- 
form of the spread of individual moment, S (Q) is the 
structure factor, given by the fourier transform of the 
distribution of these localized moments in the crystal 
and M is the magnitude of the localized moment. For 
the current loops, the magnitude of the moments is 
of order 0.1/xg, while for the spins it is 0.02/^. The 
structure factor for the current loops is proportional to 
cos (wzL) sin (2ttxq (H ± K)), where z ~ 0.29 is the ratio 
of the interlayer spacing to the lattice constant c and xq 
is the position of the centroid of the triangular plaquette 
given by (x , x ). The cosine factor reflects the fact that 
the orbital moments are identical in the bilayers while the 
sine factor arises due to the antifcrromagnetic orientation 
between the two triangles in the unit cell. For the spins, 
the in-plane structure factor cannot be completely deter- 
mined without knowing the precise relative angles of the 
spins in a unit-cell. But the fact that the spins are oppo- 
sitely oriented in the two layers implies that the structure 
factor can be as large as sin (irzL). Thus at [011] the spin 
structure factor can be as large as 1.8 times the moment 
structure factor. 

For current loops, the effective moment generated is 
spread over the area of the triangular plaquette which is 
a 2 /8, where, a is the lattice constant. For spins on copper 
and oxygen atoms, the moment is distributed over the 
d x 2_ y 2 and p x y orbitals respectively. Since the atomic 



orbitals are more localized their fourier transforms are 
weaker functions of Q as compared to the orbital mo- 
ments. To estimate the form factors we model the time 
reversal violating state with current wires along the x, y 
and —x — y directions with thickness 5. The fourier trans- 
form of this pattern of currents is expressed in terms of 
a combinations of form and structure factors. Then the 
form factor is 2exp (— TT 2 S 2 /a 2 ) /n where a is the lattice 
constant. We have assume a Gaussian profile for the cur- 
rent in the wires. Since the width of the current wire is 
related to the overlap of the copper and oxygen orbitals, 
we take it to be of order lA. Thus the form factor for 
the current loop is ~ 0.3. The form factor for the spins 
is ~ 0.9 implying that the net geometric factor for spins 
is ~ 6 times larger than those for the current loops. Give 
the estimate for the magnitude of the spin and orbital 
moments, the rough estimate for the geometric factors 
implies that indeed the resulting neutron scattering in- 
tensities due to the two orderings will be of the same 
order of magnitude. 

In conclusion, we find that given the orbital order, 
an inplane spin-order is mandated. With reasonable 
assumptions about the relative form factors for loop- 
currents and spin-moments and the calculated magni- 
tude of the ordered spin-moment, the observed polarized 
diffraction can be understood. A detailed test awaits ex- 
perimental refinements. Given the spin structure shown 
in fig. 3, DM interactions can induce a tilting of the spins 
leading to a small ferromagnetic moment. This is also 
under further investigation. 
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